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The effects of initiabstate Coulomb interactions in the charge-symmetry-breaking reaction 
dd — + QfTr" are investigated within a previously published formalism. This is a leading order ef- 
fect in which the Coulomb interaction between the two initial state protons leads to the breakup of 
the two deuterons into a continuum state that is well connected to the final ayr" state by the strong 
emission of a pion. As a first step, we use a simplified set of d and a wave functions and a plane- wave 
approximation for the initial dd state. This Coulomb mechanism, by itself, yields cross sections that 
are much larger than the experimental ones, and which are comparable in size to the contributions 
from other mechanisms. Inclusion of this mechanism is therefore necessary in a realistic calculation. 
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I. INTRODUCTION 

The concepts of charge independence and charge sym- 
metry provide powerful tools in organizing the multi- 
plet structure of systems of hadrons and nuclei. These 
symmetries are not perfect; diverse small but interest- 
ing violations have been discovered; see the reviews of 
Refs. [H, Our concern here is with charge symmetry 
breaking. 

Hadronic states can be regarded as approximately 
charge symmetric, i.e. invariant under a rotation by 
180° around the 2-axis in isospin space. Charge sym- 
metry (CS) is a subset of the general isospin symme- 
try, charge independence (CI), which requires invari- 
ance under any rotation in isospin space. In Quantum 
Chromodynamics (QCD), CS requires that the dynamics 
are unchanged under the exchange of the up and down 
quarks [l]. In the language of hadrons, this symmetry 
translates into e.g. the invariance of the strong interac- 
tion under the exchange of protons and neutrons. How- 
ever, since the up and down quarks do have different 
masses (m„ ^ rrij) [1, the QCD Lagrangian is not 
charge symmetric and neither are the strong interactions 
of hadrons. This symmetry violation is called charge 
symmetry breaking (CSB). The different electromagnetic 
interactions of the up and down quarks also contribute 
to CSB. 

Observing the effects of CSB interactions therefore pro- 
vides a probe of m„ and m^, which are fundamental, 
but poorly known, parameters of the standard model. 
For example, the light quark mass difference causes the 
neutron to be heavier than the proton. If this were not 
the case, our universe would be very different, as a con- 
sequence of the dependence of Big-Bang nucleosynthe- 
sis on the relative abundances of protons and neutrons. 
Experimental evidence for CSB has been demonstrated 



'Electronic address: 
t Electronic address: 



talahdeQu. washington.edu] 



miller@phys.washington.edu| 



many times, see e.g. Refs. [H, [3]. Two exciting recent 
observations of CSB in experiments involving the pro- 
duction of neutral pions have stimulated current inter- 
est: Many years of effort led to the observation of CSB 
m np —> d7r° at TRIUMF. After a careful treatment 
of systematic errors, the CSB forward-backward asym- 
metry of the differential cross section was found to be 
= (17.2 ± 8 ± 5.5) X 10"'' 01, where the former er- 
ror is statistical and the latter systematical. In addition, 
the final experiment at the lUCF Cooler ring has re- 
ported a very convincing dd — > an'^ signal near threshold 
(ct = 12.7 ±2.2 pb at = 228.5 MeV and 15.1 ±3.1 pb 
at 231.8 MeV) Q. These data are consistent with the 
pion being produced in an s-wave, as expected from the 
proximity of the threshold (T^ = 225.6 MeV). Studies of 
the dd QfTT^ reaction thus present exciting new oppor- 
tunities for developing the understanding of CSB. 

The reaction dd — s- an'^ obviously violates isospin con- 
servation, but more specifically, it violates charge sym- 
metry since the deuterons and the a-particle are self- 
conjugate under the charge-symmetry operator, with a 
positive eigenvalue, while the neutral pion wave func- 
tion changes sign. This reaction could thus not occur if 
charge symmetry were conserved, and its cross section is 
proportional to the square of the CSB amplitude. This 
phenomenon is unique, because all other observations of 
CSB involve interferences with charge-symmetric ampli- 
tudes. 

Due to the recent availability of high-quality experi- 
mental data on CSB, a theoretical interpretation using 
fundamental CSB mechanisms is called for. At momenta 
comparable to the pion mass, Q m^, QCD and its 
symmetries (and in particular CSB) can be described 
by a hadronic effective field theory (EFT), called Chi- 
ral Perturbation Theory (^PT), for extensive reviews see 
Refs. [7, 8, 9]. This EFT has been extended to pion pro- 
duction in Refs. [13, UH [H, [H, 1131 where typical mo- 
menta are Q ~ ^m^M, with M as the nucleon mass 
(see also Ref. [lB| where pion production was studied ne- 
glecting this large momentum in power counting). The 
EFT formalism provides specific CSB effects in addition 



FIG. 1; Diagrams relevant for the inclusion of the Coulomb initial-state interaction. The wiggly lines represent the Coulomb 
interaction between the two protons of the initial state. The dashed lines represent the emitted pion. For each Coulomb 
interaction, anyone of the four nucleons may emit the pion. 



to the nucleon mass difference. In particular, there are 
two pion-nucleon seagull interactions related by chiral 
symmetry to the quark-mass and electromagnetic contri- 
butions to the nucleon mass difference 3 13 • 

In previous work [l^, [l^ , the cross section for the re- 
action dd —> aTT^ was computed near threshold by chiral 
EFT techniques, using a chiral power counting scheme 
to assess the expected importance of different interac- 
tion terms. The first paper fis'l used a plane wave ap- 
proximation along with Gaussian bound-state wave func- 
tions. These initial calculations yielded computed cross 
sections that are a factor of ~ 2 larger than the measured 
ones. The effects of initial-state interactions and realistic 
bound-state wave functions were included later [l9j , with 
resulting cross sections of the order of several hundred 
pb or more. It is thus clear that more work is needed 
to understand the order of magnitude of the measured 
cross section, such as a treatment of the effects of loop 
diagrams along the lines of Ref. 20] . 

We consider here one specific contribution that has pre- 
viously been neglected - the influence of the Coulomb 
interaction in the initial state. The formalism employed 
is similar to that of Ref. The main objective is to 

assess the relevance of this mechanism, so we use sim- 
ple bound-state wave functions and neglect the effects 
of strong interactions in the initial state. It is worth- 
while to explain some basic features of the calculation 
that result from invariance principles: Spin, isospin, and 
symmetry requirements restrict the partial waves allowed 
for the dd — > air^ reaction. In the spectroscopic no- 
tation ^^~^^Ljl, where S,L and J are, respectively, the 
spin, orbital, and total angular momenta of the dd state, 
and / denotes the pion angular momentum, the lowest 
allowed partial waves are '^PqS and ^D^p. Hence, produc- 
tion of an s-wave pion requires that the initial deuterons 
be in a relative P-wave, with spins coupled to a spin-1 
state, forming together a state with zero total angular 
momentum. On the other hand, a p-wave pion is pro- 
duced only when the deuterons are in a relative D-wave, 
with spins maximally aligned to spin 2, requiring either 
a coupling with AL — AS = 2 or _D-states of d or a. 
Interferences between s- and p-waves disappear for any 
unpolarized observable. We shall therefore be concerned 



with the production of an s-wave pion. In the mechanism 
proposed here, the Coulomb interaction between the two 
protons of different deuterons converts the initial relative 
P-wave state into an S-wave state. Parity conservation 
then requires that one of the resulting pairs of nucleons 
be in a p-wave with orbital angular momentum unity. 
The strong pseudovector pion production operator then 
converts this p-wave state into an s-wave state. 

This paper is organized in the following manner: The 
relevant Coulomb mechanism is described in Sect. [TTl 
which also explains how the formalism of Ref. fl§\ is to 
be employed. We use the simple Gaussian bound-state 
wave functions of Ref. , but extend the calculation by 
also considering Hulthen wave functions for the deuteron. 
The detailed evaluation and numerical results are given 
in Sect. IIIII For comparison, the effects of Coulomb in- 
teractions in the final state are considered in Sec. IIVI 
These have been calculated in Ref. [l^ and were found 
to be very small. Finally, Sect. |V] assesses our results and 
discusses how these effects can be included in a realistic 
calculation that incorporates the strong interactions in 
the initial state. 



II. COULOMB MECHANISM 

In the present study, CSB arises from the initial-state 
Coulomb interaction between the two deuterons, fol- 
lowed by strong pion emission, as shown in Fig. [T] We 
note that the electromagnetic contributions can be or- 
dered [l^ relative to each other in the same fashion as 
the effects of strong CSB. In this case, the leading or- 
der (LO) term considered here is of O [a^^M / {in f^p)] , 
with M as the nucleon mass. This term is roughly of 
the same size as the LO strong CSB term which is of 
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^m^/(M4,r/3p)J. 
The CSB pion production operator Oq is given by 



O. 
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O, {E-H„ 
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where is the Coulomb interaction between the two 
protons in the initial state, which acts to form a four- 
body continuum state that propagates according to 



3 



{E — Hq + ie) ^. It is convenient to write as a sum 
of pair- wise operators: 

Vc -- 



i.k 
k^c > 



(2) 



where the Qj k nucleon charge operators. The oper- 
ator Hq is the sum of the kinetic energies of each of the 
four nucleons. The strong pion production operator is 
denoted by ^^'^ is given by 



fe, 



2f, I ^rl^ 
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where the fcjjfcj are nucleon momenta before and after 
the pion emission, respectively. The p-wave term with 

= — can be ignored in the threshold regime consid- 
ered here. The factor w is also replaced by the pion mass 
[I = 134.974 MeV. 

The present analysis uses a plane- wave approximation 
and the simplest possible d and a bound-state wave func- 
tions, those of a Gaussian form. Assuming spatially sym- 
metric bound-state wave functions, the invariant ampli- 
tude is given by 



with 



M 



\A) 
\DD) 



cPr(fipid% {A\0\DD), 



V^*d(Pi)*d(P2) \dd), 



(4) 



(5) 
(6) 



where and <i>^ are the spatial parts of the a-particle 
and deuteron bound-state wave functions, and s = AE^ 
is the total cm. energy squared. The ket vectors \a) and 
\dd) contain the fully anti-symmetrized spin and isospin 
wave functions. Because of symmetry requirements, the 
plane- wave dd scattering wave function is included in \dd) 
as given by Eqs. and ITTl) below. The invariant am- 
plitude can then be written as 

M = y^2E^s J d^rd\d%^i{r,p,,p^) 

{a\0\dd) <i>MMP2), (7) 

The matrix element {a\0\dd) contains all the spin- isospin 
couplings of the nucleons and the pion production opera- 
tor O. The wave functions are expressed in terms of the 
(2-1-2) Jacobian coordinates 

R = ^(ri + r4) (= in cm.), 



Pi = ri-r^, 



P2 



3 ' 4' 



(8) 



with the corresponding momenta 
K 



fc]^ + ^2 + '^3 + (= in c.m.), 
1 

2^ 
1 

2* 
1 

2^ 
1 

2* 



k — -,(^1+^2 ^3 ^^4) 

j(Pi-P2) (=pinc.m.), 
-{ki — fej), 
t9 = ^(^3 — k^), 



(9) 



defined so that ^i'^i = K ■ R + k r + k^ - pi + - P2- 
The Jacobians are equal to unity in both representations. 

The ground-state wave functions are represented by 
Gaussian functions, and these may be explicitly ex- 
pressed in the above coordinates using X)i<j(^i ~ ''"j)^ = 
4r^ + 2pl + 2pl, yielding 

*a(^,Pl>P2) = 
^diP^) = 



^9/4q,9/2 
1 

7r3/4/?3/2 



exp 



{2r^+pi + pl) 



(10) 



exp 



PL 

2/32 



z = l,2 (11) 



where the parameters a = 2.77 fm and /? = 3.189 fm 
have been fixed using the measured a and d rms point 
radii (r2)i/2 = 1.47 and {rjy^^ = 1.953 fm [21^. We 
shall work in momentum space and therefore record the 
corresponding wave functions 



exp 
87r9/4 ' 



^'a(fc,Ki,/«2) 

^dii^t) = exp ( - 

TT 



i^l,2 



3/4 



(12) 



(13) 



In order to study the sensitivity of our results to the 
choice of wave functions, we also use a deuteron wave 
function of the Hulthen form: 



c&f (kJ = TVf 



1 
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62 



i = 1,2 



■\/ ab{a + b) 
7r(a — b) 



(14) 



where the parameters are given by a = 0.23161 fm ^ and 
b= 1.3802 fm-i [13]. 

Since we have assumed that the orbital parts of the 
wave functions are symmetric under the exchange of any 
pair of nucleons, we may define the initial- and final-state 
spin-isospin wave functions as 
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1 



V2 



{((l,2)i,(3,4)i)J[l,2]o,[3,4yo 



4 



((l,2)o,(3,4)o)o[[l,2]i,[3,4]Jo}. 



(I«l> + I«2>) 



\dd) 



Ml2'^34) 



1 

^ (1 - P23 - A4) Ml2'^34> , 

((1,2)1, (3,4)i)s[[l,2]o,[3,4]o]o 



(15) 
(16) 

(17) 



where and [«, jlj^ are the spin and isospm Clebsch- 

Gordan couplings, with magnetic quantum numbers sup- 
pressed, for nucleons, or nucleon pairs, i and j couphng 
to spin s and isospin T, respectively. We shall refer to 
the first term of Eq. (jlSp as the "dd" component of the 
a because the pairs (12) and (34) each have the spin and 
isospin of the deuteron. In the above equations, P^j is 
the permutation operator of the indicated nucleons. The 
symmetry requirements for the exchange of the deuterons 
are represented by the (orbital angular momentum de- 
pendent) combination of plane waves in Eq. (|17p . with 
p as the relative momentum of the deuterons. Even 
though the expression for the a state superficially sin- 
gles out a (12)-|-(34) configuration, it is indeed fully anti- 
symmetric in all indices. This particular form is used 
because it closely matches the form of the initial-state 
wave functions, thereby simplifying the evaluation of the 
spin-isospin summations in the matrix element. In prac- 
tice, the dd wave function can be simplified to 

\dd) = V6((l,2)i,(3,4)i)5[[l,2]o,[3,4]o]„ e^P-^, (18) 



since each of the three terms in Eq. pB]) gives an identical 
contribution to the matrix element. 

The expressions (flS)) through provide insight that 
simplifies the calculation: The Coulomb interaction has 
no spin operator, so the initial state is connected only 
with the "dd" component of the alpha particle, which 
means that only the Class III [l| part of the Coulomb op- 
erator (r^ j + '''z.j) contributes. One of these operators 
finds another from within the pion production operator 
and is squared to ± unity. As a result, the matrix ele- 



ment turns out to be proportional to the spin operators 
of the (1,2) and (3,4) systems. In the normalization used 
here, the spin-averaged cross section (for s-wave pions) 
is given by 



16ns \p\ 9^' ' ' 

pol. 



(19) 



where the summation is over the deuteron polarizations. 



III. EVALUATION 

The analysis is most conveniently performed in mo- 
mentum space. Combination of Eq. ([T]) with Eq. ^ 
yields, upon Fourier transformation, the Coulomb con- 
tribution of present interest, A4(j '■ 

X (a,|$.(fc,K„K,)^ ^^Z:^ — 
^=i'4S — 



X {k,K^,K2\Vc\dd,p) 



ki 



2k\ 



2k,2 



(20) 



fe has been 



where the relation X]j=i 4 '"j 
applied in the free propagator, and \dd^p) represents the 
initial dd (relative plane wave) state of Eq. ((T8)) . includ- 
ing the internal spatial, spin, and isospin degrees of free- 
dom. Further, (fc, Kj^, /«2|V'p|d(i,p) denotes the momen- 
tum space representation of the state formed by the ac- 
tion of Vf-, on the initial state. The round bracket nota- 
tion used here signifies that only the spatial degrees of 
freedom are included. 



A. Reduction to Quadrature 

The first step in the calculation is the simplification 
of the pion production operator in Eq. (|20p . Define the 
operator X according to 
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^ Tza) 
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(22) 
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It is instructive to consider the first term of Eq. ([22]) : It should be noted that the operators Vq^, Vq^ do not flip the 
spin of their deuteron. Also, the initial-state deuteron (1,2) is connected to the deuteron-like (1,2) component of 
the a. Thus the terms with i,t^ 2 '^^^ be dropped, and the initial-state deuteron (3,4) is similarly connected to 
the deuteron-like (3, 4) component of the a. We are required to have CSB, so only the terms proportional to 3 are 
relevant. For the (3,4) "deuteron" of the final state we have 3 = — 4. Thus the first term of Eq. ((22)) simplifies 

to j{cr^ ■ — cr^ ■ k^)(v'^^ + f^'^). Similar manipulation of the remaining terms in Eq. (|22p leads to the result 



X 



"■3 -^3 



^4. ' ^4 



3.2 4,1 



3,1 



3,2 



(23) 



Next define the spin operators S-^ = \{cti + cTr^iS^ — \{<y-j, + cr4),Sj^ — ^(f 4 — cr2),S2 = \{p''i ~ "'4), such that 
each of the cr^ is a linear combination of the and Sj. Only the terms proportional to connect the initial state 
to the "dd" component of the a. Thus one finds 



X 



S2 ■ K2 



3,1 



Si ■ Ki 



3,1 



3,2 



(24) 



Further, it is permissible to interchange indices 3 and 4 in the spatial wave functions multiplying the first term 
of Eq. (j24p , and similarly to interchange 1 and 2 in those multiplying the second term. The final form of the operator 
X is thus 



X 



S2 • 



3,1 , 3,2 



Si ■ Ki 



3,1 



4,1 



The next task is to compute the momentum space matrix element of the operator X . This is given by 



{k,Ki,K2\X\dd,p) = S2-i^2 {k,i^i,i^2\{^c^ +-^c'^)\dd,p) + Si • Ki (fc, /«4, KjI (t^, 
with the spatial matrix elements 



3,1 



\dd,p) 



{k,ni,n2\v^'\dd,p) = M^PArf^^ ^^.^,p,^.^,.p,-.k.r ^p.r a>^(p^).,^(p^), 

J (27r)y/^ \r,-rf^\ 



(25) 
(26) 

(27) 



which are thus found to be products of the momentum-space Coulomb interaction with deuteron wave functions 
evaluated at shifted values of the momentum. In particular, we define v = k — p and obtain 

{k,Ki,K2\v^c'W,p) = (IJP^ (/^l - f ) 5d (/^S + I) , 

{k,Ki,K2\vl^\dd,p) = (^^J^s/T^, (/^l + I) (/^2 + I) , 

{k,K,i,K2\v'^'\dd,p) = ^ipl^ (/.^ _ I) (;.2 - I) • (28) 

Insertion of these results into Eq. (pS)) finally gives 

87ra, 



{k,Ki,K.2\X\dd,p) 



(27r)3/2^;2 



(^2 • ^^2 - "^1 • *d «1 + 77 «2 + 77 



(29) 



Examination of Eq. (j20p reveals that p is the only momentum remaining after the integrals have been performed. 
Thus the terms with Ki and k,2 in Eq. (|29p both end up being proportional to p, which we may take as the z— axis. 
Furthermore, as the integrands of both terms are identical, the whole operator must be proportional to (5*2^ — S'l^) 
and the integrand to ^^'t'^^ . We thus need to consider the spin matrix element 



(lAfi, IM2IO, 0)(lMi, 1M2|(52, - ^iJllMi, lAf2) - (IMi, 1M2|0,0)(M2 - M^) 

= 2M2(lMi,lM2|0,0), 



(30) 



where {IMi, 1M2|0,0) is the Clebsch-Gordan coefficient that couples the spins in the "dd" component to zero. Armed 
with this knowledge, we may now use Eqs. (HH) and (|29I30|) in the matrix element (HO]) to obtain 



Mr 



U J (2^)3/2 



X M2(lAf4,lAf2|0,0) 



(31) 



d^kd^Kid^K2 



{1^1 + 1^2 ) -P 



[2ME - {2kI + 2kI + fc2) + It] ' ^^{k, Ki,K2) *d (^i + I) *<i ('^2 + I) 
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The above equation is our main result, and it allows for the use of general radial wave functions. However, because 
of the zeros in the energy denominator, it may not be well suited to evaluation using Monte-Carlo techniques. 
Nevertheless, if certain simple wave functions are used, Eq. ([3T|l may be simplified further. In the next subsection, 
this will be performed for wave functions of the Gaussian and Hulthcn types. 



B. Gaussian and Hulthen Deuteron Wave Functions 



If the Gaussian wave functions of Eqs. (fT2|) and (fT3| are used, the expression (j3ip becomes 

fig A Sira^ 



L J (2^)3/2 

3 + K^)-P 



X M2(lAfi,lM2|0,0) X N^Nfi x /„, 



d^kd^K^d^K^ ' [2ME - {2kI + 2kI + k^) + ie\ 



X exp 



2 , r,._2\ 



{k^ + 2kI +2k^) 















exp 




'2) _ 


exp 




'2) _ 



(32) 



(33) 



where we recall the definition v — k — p. The factors in the denominator of Eq. ([33]) may be rewritten in terms of 
Gaussians, giving 

POO 

v-^ = dj exp [-7(1;^ + ei)] , (34) 
Jo 

[2ME - {2kI + 2kI + k^) + it] = -i dv exp [ iv (2ME - (2^1 + 2k\ + k^) + ie)] , (35) 

Jo 

where the regulator e assures that the integral over v converges, and is included to handle the point k — p. The 
use of the above identities leads to an 11 dimensional integral, of which 9 dimensions involve products of Gaussian 
functions, such that these integrals may be computed analytically by successive completion of squares in the exponents. 
This procedure yields a two-dimensional integral over v and 7, which is then computed numerically. In this way, using 
the definitions 



the integral /„ may be re-written as 



+ 1, 



I 2 I 



21' 
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Y 



(36) 



Ig = —i I dv d'y I d vd nd I up exp 



- + -2 



X exp 



2 02 



7('y'+ei) 



2 4 

where the completion of squares is facilitated by the definitions 



exp {iv (2ME - 4i^ - /t^ - (t> + pf + it)\ , 



= -I- 8w, R} 



^+/32 + 4zi., Rl^\ 



P 



P 



2M 



(37) 



(38) 



For all successive equations, we will define Ij — —ilpln^^^Iij. At this point, the regulators e and ei may be safely 
dropped. By combination of the above results and definitions, we find for the case of Gaussian deuteron wave functions 



- [- — ) J^2 



X M2(lMi,lM2|0,0) X N^Ni x {~-i\p\n^/^) x 1^^, 



, , J 2 

'^'^^^;r8^ ^""Pv p 



exp 



ivp a 



(39) 



If the Hulthen wave function, given in Eq. (|14p . is used for the deuteron, the product of deuteron wave functions 
in Eq. pip may be re-written as 



dr]idrj2 \ exp 



X < exp 



-V2 I '^2 + 2 



(40) 
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FIG. 2: Diagrams relevant for the inclusion of the Coulomb final-state interaction. The wiggly lines represent the Coulomb 
interaction between the two protons in the final state. The dashed lines represent the emitted pion. 



after which the calculation proceeds, as before, through successive completion of squares, but in the Hulthen case 
we are left with a four-dimensional integral suitable for numerical evaluation. The end result is that the integral Iig 
of Eq. ([59)1 should be replaced by 7^^, which is given by 



J-j I 



8Rl)^ Rl 



exp ( ) fab{Vl,V2), (41) 



where the definitions 



Vi2~ ] ' -K, - — -t- -I- ?72, Rk--^+Vi2^ -Ku - -K^ +7- — , (42) 



fabiVi,V2) = e-'''(''i+''2) +e-'''(''i+'"2) -2e-"'''le-'''''^ (43) 

are used. These expressions represent the complete amplitude. In order to obtain the cross section, it is necessary to 
evaluate the spin sum 

M|(lMi,lAf2|0,0)' = I (44) 



and insert everything into Eq. (|19p . yielding 

. \2 |„||„ U2 ^9/:j6 

(45) 

(46) 

where dg and cr^ again denote the expressions relevant for the Gaussian and Hulthen deuteron wave functions, 
respectively. It should also be noted that the above expressions remain valid for the Coulomb interaction in the ^He 
bound state, which is considered in Sect. IIVI 



- "^HTTJ 288 I'l^l ' 

' ~ "[17) 288 ^ (a - 5)4 l^i'^l 



IV. COULOMB IN THE *HE BOUND STATE 



A complete assessment of all Coulomb effects should include a treatment of the Coulomb interactions in both the 
initial and final states. Our focus here is on the initial-state effects, as these have not been considered up to now. 
However, it is also worthwhile to compute the effects of the Coulomb interactions in the final state within the present 
framework. We recall that in this framework, the strong interaction between the initial-state deuterons is neglected, 
and simple bound-state wave functions are used. If the effects of Coulomb interactions in the ^He bound state are 
included, the reaction dd — > a7r° can proceed via strong pion production, which is here assumed to be initiated by 
the one-body operator. The relevant CBS pion production operator is then given by 



(47) 



8 



where ~ 28.3 MeV is the ^He binding energy. It is instructive to define the operator Y according to 

i,j<k=l,4 



(48) 



fl + r. 



(l + ^.,3) 



3,2 



(1- 


f T.,2) 




2 


(1- 






2 


(1- 


^ ^2,4) 




2 


(1- 


^ ^2,4) 



3,2 
c 


(t^,3'^3 


^3 


+ 'rz,4'^4 


^4) 


4,2 
c 


(^^,3*^3 


fc3 


+ ^2,4*^4 


^4) 


4,1 

C 


(Tz.iC^I 


fcl 


+ ^;^,2'^2 


^2) 


4,2 
C 


(^^.I'^l 




+ ^3,2*^2 


^2) 



(49) 



In the above equation, the one-body operator can either maintain the spin-parity-isospin quantum numbers of a single 
deuteron, or produce a single two-nucleon state with 5* = 0. As the operators v^^,Vq'^ do not flip the spin of their 
two-nucleon system, the quantum numbers of the dd state must be maintained. Manipulations similar to those of 
Sect, mil lead to the simplification 



Y 



3,1 , 3,2 4.1 4,2 
V -^c -^c 



(T3 • ^3 — (T4 • ^4 



3,1 , 4,1 3,2 4,2 



and finally to 



Y 



3,1 , 3,2 



3,1 , 4,1 



Si ■ 



, (50) 



(51) 



which is analogous to that of Eq. (j25p for the initial-state Coulomb interaction. Computation of the momentum space 
matrix element then leads to the result 



{k,Ki,K2\Y\dd,p) = 



(27r)3/2t;2 



82 - {^2 + -] - Si - {ki + - 



(52) 



Since p is again the only momentum remaining after integration, the above matrix element may be treated along the 
same lines as Eq. ([29)) . Thus it is again possible to extract a factor (5*22 — •S'l^), giving finally 



M-c = (-^) 



(27r)3/2 



l; = / d-^kd-'K^d-'K. 



(Ki + K2+v) -p 



a 

2^-■-l 



X exp 



2 



[-2MeB - {2kI + 2kI + k^)] 

( V 

2 V 2 ^ 2 











exp 


[4 (".-1)1 


exp 


2 



(53) 



where the Gaussian wave functions of Eqs. (fT2|l and (fT3|) have been employed. We proceed by writing the last two 
factors of Eq. ([33]) in terms of Gaussians, and note that the only difference with the general procedure of the previous 
section is that we may use 

1*00 

[-2M€b - {2kI + 2kI + k^)] = - dv exp [-v (2Af -I- 2k\ + 2kI + k^)] , (54) 

"'0 

upon which Ig becomes an 11-dimensional integral, of which 9 can again be computed analytically by completion of 
squares in the exponents. Analogously to the previous section, we employ the notation 



2 4 4Rf 



(55) 



along with the definition = |p|7r^/2/jF. ^j^g matrix element in Eq. ([53| then becomes 



dv d'y 



- 



/. ; (27r)3/2 
exp 



X M2(lMi,lM2|0,0) X N^Nj x \p\tt'>'^ x /f^, 



54 



2i?fi?5 



1 - 



8i?2/ 



exp {—2v Meg) . 



(56) 
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The evaluation of the final-state Coulomb mechanism of the preceding subsection can also be implemented using 
the Hulthen wave functions for the deuteron, given in Eq. (dH). The net result is that the integral if^ should be 
replaced by with 

= Jo '"'^ '''''''' ^iwm 

where the definitions 

C12 - ' -^i - Y ^1 -^K~"^ + ''i2' ^v-^K+l~^^ 

are used. The cross sections can then be computed using the expressions given in the previous section on the initial- 
state Coulomb interaction. , 



exp(-2z/MeB) /afc(f7i, 



(57) 



V. NUMERICAL RESULTS AND DISCUSSION 

The model parameters used in the present calculations 
are given in Table U and the calculated cross-sections at 
= 228.5 MeV and = 231.8 MeV are summarized in 
Table |TT1 such that and ct^ denote the results for the 
initial-state Coulomb interaction, obtained with Gaus- 
sian and Hulthen deuteron wave functions. If Gaussian 
wave functions are used throughout, the results are 59 pb 
and 75 pb at the two energies considered. If Hulthen 
wave functions are used instead for the deuterons, these 
results increase to 87 pb and 111 pb. In either case, the 
effects of the initial-state Coulomb interaction are signif- 
icant, as the experimental values are 12.7 and 15.1 pb, 
respectively. It should also be noted that these differ- 
ences are much smaller than those encountered between 
the CD-Bonn and Argonne VIS potentials in Ref. p^ . 
The first toy-model calculations yielded nominal values 
of 23 pb and 30 pb. The present mechanism is there- 
fore clearly large enough to warrant inclusion in a fully 
realistic calculation. 

TABLE L Summary of parameters used in the calculation. 
The values of a and /3, which appear in the expressions for 
the Gaussian bound-state wave functions, are from Ref. [2H . 
whereas a and b are relevant for the Hulthen deuteron, and 
have been taken from Ref. [2^ . 



a [fm] 


2.770 


13 [fm] 


3.189 


a [fm"^] 


0.23161 


b [fm-^] 


1.3802 


Mes[fm"^] 


0.68 


U [MeV] 


92.4 


9a 


1.26 




137.04 



The results for the Coulomb interaction in the final 
state are denoted and cr^, and are also given in Ta- 
ble [Til If Gaussian wave functions are used throughout. 



TABLE II: Summary of momenta, calculated cross-sections 
and integrals for the matrix elements. The quantities with the 
subscripts g and h correspond to the Gaussian and Hulthen 
deuteron wave functions, respectively. The energies and mo- 
menta correspond to those of Ref. [6|]. 





T<i = 228.5 MeV 


Ta = 231.8 MeV 


p [MeV] 


462.913 


466.924 


P. [MeV] 


19.372 


28.266 


[MeV] 


3727.430 


3727.487 


ha [fm-5] 


(1.3638, 1.3919) x 10"^ 


(1.2440, 1.3213) X 10"^ 


Kg [fm-^] 


2.3043 X 10"^ 


2.1713 X 10"'' 


hu [fm-^] 


(2.7011,2.1617) X 10'=* 


(2.4919,2.0821) x 10"^ 


I(h [fm-^] 


1.1465 X 10~^ 


1.0957 x 10"^ 


O-g [pb] 


58.95 


75.25 


cffe [pb] 


85.77 


111.2 


[pb] 


0.824 


1.077 


[pb] 


9.419 


12.66 



the results are much smaller, about 1 pb, which repre- 
sents '--^ 1% of those found for the Coulomb interaction 
in the initial state. The use of Hulthen deuteron wave 
functions in the initial state is found to enhance the ef- 
fects of the final-state Coulomb interaction. However, 
they are still relatively small, about 10% of those of the 
initial-state Coulomb interaction. 

The principal result of this study is the manifest need 
to incorporate the effects of Coulomb interactions in the 
initial state into the realistic calculation that includes 
strong initial-state interactions. Recent progress in the 
treatment of Coulomb interactions in few-body scattering 
calculations [23] should eventually allow such computa- 
tions to be performed. 
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